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^ '■ 1 General Introduction 



The Journal Club Talk "Dark Matter in Galaxies" was delivered within the "Dark Matter Awareness 
■ Week (1-8 December 2010)" at 140 institutes in 46 countries and was followed by 4200 people. More in- 
' formation on DMAW can be found at its website: www.sissa.it/ap/dmg/index.html and a short movie 
• • ' that documents this worldwide initiative at: http : //www. youtube . com/watch?v=A0Bit8a-lFw, 
. ^ ' For this DMAW 2010 talk there was a standard reference presentation with explaining notes 

^ , prepared by a coordinated pool of leading scientists in the field, which the various speakers adapted 
I and delivered in their Institutes so as to suit local circumstances, audiences and their own personal 
preferences. Importantly, the speakers subsequently provided us with comments on the reference 
material and on the initiative itself (see Section 1.3 to download them). 

After the event, in response to participants' feedback, we improved and upgraded the guideline 
material to a "Presentation Review" on the subject "Dark Matter in Galaxies", for which we provide 
a .pptx file and the related Explanatory Notes. This is described below together with suggestions of 
who might be interested in using it and how to do it. 



1.1 What is the Presentation Review + Explanatory Notes? 

This is an innovative scientific product, introduced for fostering scientific debate on Cosmology and 
for rapidly diffusing, among the scientific community, results in a field where, we believe, knowledge 
is rapidly growing. Supported by the large number of speakers delivering the DMAW 2010 talk, by 
their reactions and by those of their numerous audiences, we decided to address the whole community 
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of cosmologists and astrophysicists. We prepared this Presentation Review, convinced that there is a 
generahzed request of knowledge about the DM mystery at galactic scale. 

The Presentation is a Review Talk at post-graduate level, (2 x 45) minutes long, on the specific 
subject: "The distribution of Dark (and Luminous) Matter in Galaxies", although other cosmological 
issues (e.g. galaxy formation) are treated, but from the phenomenological perspective of the galaxy 
mass distribution. 

The Presentation Review is not a PowerPoint translation of an already existing Review Paper, 
but it has been prepared directly as a Presentation on the above described subject. The Explanatory 
Notes are not the written version of the Presentation, but concepts, formulas and references that help 
with understanding the slides and preparing an eventual speaker to comment on them to an audience. 

In detail, their scope is to enable the user who reads them to follow the Presentation easily and 
to get the greatest benefit from it. The user who, instead, chooses to carefully study them, would be 
able herself/himself to deliver the Review Talk. For this reason, there is some additional material in 
the Notes (arguments, references, formulas) for the benefit of a prospective speaker or any user who 
wants to deepen the knowledge on the subject. 

The presentation runs also two movies in prefixed slides. These movies are especially devised 
for the Presentation and concern the Cosmological Simulations and The Universal Rotation Curve of 
galaxies. 

It may be obvious, however, it could still be worth to stress that the Explanatory Notes without 
the .ppt Presentation are worthless. On the other hand, the Presentation Review file without the 
Notes is difficult to "use". 

An other clarification maybe necessary: we are envisaging this Review as being one section of a 
virtual course on "Physical Cosmology" organized as follows: Big Bang Cosmology (3 x 45'), Structure 
formation (4 x 45'), Dark Ages and objects at high redshifts (3 x 45'), The distribution of DM in 
Groups, Clusters and larger scales (2 x 45'), The distribution of DM in Galaxies (2 x 45'), The 
physics of baryons in forming the luminous structure of Universe (4 x 45'), Candidates for the DM 
particle and their detection (3 x 45'), Alternatives to Dark Matter (1 x 45'). We have not written 
out this programme to anticipate a future project, but to frame this Review Talk in its Cosmological 
Context. More specifically, our Review aims to bring the user up to a good level of knowledge of its 
core subject: the Dark and Luminous matter distribution in galaxies. For other topics, which are 
occasionally touched on by the review, the user is redirected to more specialized reviews, lecture notes 
and books. And, of course, a number of issues involving DM will not be dealt in this Review. 

Let us also stress that, contrary to the situation for written reviews, with this one it is up to the 
users, and not only the authors, to establish the relative and absolute importance of each statement, 
formula, reference and plot in the Presentation and eventually convey this to others. 

1.2 Who may this Presentation Review + EN be useful to? 

In order to avoid disappointments, it is worth to stress that this Review is neither designed as a public 
outreach talk on Dark Matter, nor as a Seminar that advertises a hot topic in Cosmology. 

The target of this Presentation Review is the "Community" of thousands of physicists and astro- 
physicists whose work already touches, in many different ways, the subject of dark matter in galaxies 
and that would like to become more aware of the phenomenology behind this subject, possibly relevant 
for their research. 

The success of DMAW 2010 has shown that this community is wide and that our dedicated efforts 
are well received. We hope that many people will be interested in downloading this Presentation + 
the Explanatory notes: 

i) People who would like to follow this innovative approach so as to learn more about this increas- 
ingly important subject in Cosmology, Extragalactic Astrophysics and Astroparticle Physics, which 
they may often meet with in their own research work. 

ii) PhD Students in any field of Astrophysics and Astroparticles 

iii) Newcomers in Cosmology from different research fields. 
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iv) People who would like to deliver a well-tested Presentation Review to colleagues or students, 
for the reasons given above. 

v) People that would like to have a reliable Review Presentation on an important subject that may 
be of help in the preparation of talks or lectures, centred either on this or on other adjacent topics. 

1.3 Downloading 

The Presentation Review, the Explanatory Notes, the two movies that run inside the Presentation 
and the Comments received for the DMAW 2010 Initiative can be downloaded at: 

|http : //www. sissa. it/ap/dmg/dmaw_presentation.html 

In detail, the names of the files of the downloadable material are: 

The Presentation Review Talk "Dark Matter in Galaxies": DMGAL2010.pptx 

Explanatory Notes: EXPLANOTES.pdf 

Movies in the Presentation: moviel.mp4, movie2.wmv 

The DMAW2010 Speakers' Reports and Comments: Comments.pdf 

1.4 How to use this Presentation Review + EN ? 

The following may be obvious but, to avoid misunderstandings, must be stressed for a correct and 
successful use of this material. The use of this innovative product requires significantly more time and 
personal effort than reading a Review Paper published in a Journal or in a Conference Proceedings 
Book. Users with profile (i-iii) above will need to invest and work a couple of half- days to read 
the Notes and to play the Presentation. For users with profile (iv) the requested effort is obviously 
greater. For users with profile (v), instead, the material we provide can save their working time and 
enable them to get improved results with less effort. Finally, let us talk about the users^, that is, those 
who will be in the audiences of the users with profile (iv). They are the best off: with 90 minutes of 
attention they will get a good knowledge of the field, as it was the case for the audiences of DMAW 
2010. 

1.5 How to refer to/cite this Presentation Review + EN ? 

Differently from the usual situation for reviews published in Journals and Books, the present new 
scientific product should not be cited in papers as the unique support reference for a result or for a 
view, also if they are its core subject. We believe that the credit to the results that this review wants 
to diffuse should go to the original papers. On the other hand, we would much like to ask people to 
advertise this Presentation Review as much as possible in papers and in other occasions (e.g. Social 
Networks etc), in terms of the aims described in Section 1.2 

1.6 Comments and criticisms 

They are welcomed, especially if they consist in modifications and additions in a form that, in case, 
allow us easily to take them into account. 

1.7 Acknowledgments 

We warmly and sincerely thank our colleagues, all the speakers, the LOC and the SOC because, 
without their willingness, suggestions and support, this product so as the whole DMAW 2010 initiative 
would not have been realized. 

cfrigerio@ufabc.edu.br; Andrea.Lapi@roma2.infn.it, salucci@sissa.it 
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2 The Review Talk: Dark Matter in Galaxies 
Abstract 

This Review Talk concerns, in a detailed way, the mass discrepancy phenomenon detected in galaxies 
that usually we account by postulating the presence of a non luminous non baryonic component. In 
the theoretical framework of Newtonian Gravity and Dark Matter Halos, we start by recalling the 
properties of the latter, as emerging from the state-of-the-art of the numerical simulations performed 
in the current ACDM scenario of cosmological structure formation and evolution. We then report 
the simple, but much-tclling, phenomenology of the distribution of dark and luminous matter in 
Spirals, Ellipticals, and dwarf Sphcroidals. It will be shown that a coherent observational framework 
emerges from reliable data of different large samples of objects. The findings come after applying 
different methods of investigation to different tracers of the gravitational field. They include RCs 
and dispersion velocities profiles fitting, X-ray emitting gas properties analysis, weak and strong lens 
signal mass decompositions, analysis of halo and baryonic mass functions. We will then highlight the 
impressive evidence that the distribution of dark and luminous matter are closely correlated and that 
have universal characteristics. Hints on how this phenomenological scenario of the mass distribution 
in galaxies, including the Milky Way and the nearby ones, has a cosmological role, are given. Finally, 
we discuss the constraints on the elusive nature of the dark matter particle that observations pose to 
its direct and indirect searches. 
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3 Explanatory Notes for the Slides 

Introduction 
Slide 2 

Dark Matter is a main actor in Cosmology. It rules the processes of formation and evolution of all 
the structures of the Universe, of which today it constitutes the great majority of the mass and an 
important fraction of the mass energy. Furthermore, it is likely made of an elementary particle, beyond 
the Standard Model, but yet not detected. However, this review will focus on just one specific aspect 
of this mysterious component, although a very crucial one: its distribution in galaxies, examined also 
in relation with that of the ordinary matter. The phcnomcnological scenario discussed in this Review 
Seminar can result a precious and unique way to understand how dark halos and their galaxies formed 
and what they are really made of and maybe approach ourselves to new laws or processes of Nature. 

Slides 4 

Galaxies, i.e. Spirals, Ellipticals, and Dwarfs spheroidals, have very likely just one type of dark halo. 
Their luminous components, instead, show a striking variety in morphology and in the values of their 
structural physical quantities. 

The range in magnitudes and central surface densities is 15 mag and 16 mag/arcsec^. 

The distribution of the luminous matter is given by a stellar disk + a spheroidal central bulge and 
an extended HI disk in spirals and by a stellar spheroid in ellipticals and dSphs. 

Slide 5 

The best physical way to introduce the "phenomenon" of Dark Matter is the following: let us define 
M(r) the mass distribution of the gravitating matter and Milr) that of all baryonic components. 
We can obtain both distributions from observations; we first realize that, as we start moving from 
a certain radius vt, function of galaxy luminosity and Hubble Type, outwards, these distributions 
increasingly do not match, c/logM/dlogr > d log M^/c/ log r. 

We are then forced to introduce a non luminous component whose mass profile Mjj(r) accounts 
for the disagreement: 

d log M(r) _ML{r)d log Ml Mnir) d log Mh 
dlogr M{r) ci log r M(r) dlogr 

This immediately shows that the phenomenon of the mass discrepancy in galaxies emerges from 
the discordant value of the radial derivative of the mass distribution, more precisely, from that of the 
circular velocity or that of the dispersion velocity, two quantities that measure it: M oc r y^(r) oc 
a2(r). 

The DM phenomenon can be investigated only if we know well the distribution of luminous matter 
and we can accurately measure the distribution of the gravitating matter. 

Slides 6 

We briefly present the aspects related to the current theory of galaxy formation that are mostly 
relevant for this core topic of this Review. 

Slides 7-8 

It is worthwhile to briefly recall the predictions for the structure of DM halos coming from ACDM, 
the current theory of cosmological structure formation which is in well agreement with a number of 
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cosmological observations. 

A fundamental prediction, emerging from N-body simulations performed in this scenario is that 
virializcd dark matter halos have a universal spherically averaged density profile, Pif(r) = pcDM{f) 
(Navarro, Prenk, & White, 1997): 

PCDMyr) = , . l—yi' 

( M 

where ps and Vg are strongly correlated (Wechsler et al., 2002, 2006), ^ 8.8 [iQif^j kpc. The 
concentration parameter c = Tf^/ Rvir '^^ found to be a weak function of mass (Klvpin et al., 2010) but 
a crucial quantity in determining the density shape at intermediate radii. 
Notice that for any halo mass distribution Mff(r): M^r = MniRvir) 

The recent Aquarius simulation, the highest resolution to date, renders a single dark matter halo 
using 4.4 billion particles, of which 1.1 billion within the virial radius. We show spherically averaged 
density (left) and circular velocity (right) profiles. Curves with different colours correspond to different 
resolution runs, and the corresponding resolution limits are highlighted with arrows. These show that 
the spherically-averaged density profiles of ACDM haloes deviate slightly but systematically from the 
form found by NFW by using a much lesser number of particles and become progressively shallower 
inwards. However, these resolution differences are quite irrelevant and the NFW circular velocity well 
represent ACDM halos. Notice that the latter is almost never flat. 

In any case, the actual profile is of very uncertain origin: departures from NFW shape are also 
seen in the halo velocity dispersion profile a and the pseudo phase-space density emerges as a power 
law of radius: p/a^ oc r^-^ e.g. (Lapi &; Cavaliere, 2009) 

Slide 9 

The movie shows the formation of one of the Aquarius halos, over nearly the full age of the Universe. 
The camera position moves slowly around the forming galactic halo, pointing towards its centre at all 
times. 



4 Dark Matter in Spirals 
Slide 11 

The stars are distributed in a thin disk with surface luminosity (Freeman, 1970): 

where Rd is the disk length scale, Iq is the central value and M/j is the disk mass. The light profile 
does not depend on galaxy luminosity and the disk length-scale Rd sets a consistent reference scale 
in all objects and the extent of the distribution of the luminous component, we take Ropt = 3.2Rd as 
the stellar disk size. 

Although the mass modeling of some individual object and certain investigations require to accu- 
rately deal with the occasional presence of non exponential stellar disks and /or with the quite common 
presence of a central bulge, the above equation well represents, for the aim of this talk, the typical 
distribution of stars in spirals 

Slide 12 

Spirals have a (HI + He) gaseous disk with a (very) slowing decreasing surface density. This disk, at 
outer radii R > Ropt, especially in the case of low luminosity objects, is the major baryonic component 
and therefore must be carefully considered in the mass modeling. An inner H2 disk is also present but 
it is negligible with respect to the stellar one (Wong Sz Blitz, 2002) 
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Slide 13 



The kinematics in spirals (and more rarely in other objects) is measured by the Doppler Shift of 
well-known emission lines of particular tracers of the gravitational potential: HI, CO and H^. 

RadioTelescospes measure the galaxy kinematics by exploiting the 21 cm emission line. Notice 
the importance of the spatial resolution: since if we want to map a galaxy with a sufficient number 
of independent data also in the inner regions, the HI emission disk size must be at least 15-20 times 
bigger than the radiotelescope resolution beam, i.e. ~ 200". There are relatively few objects in the 
Universe complying with this requirement. 

WSRT (Westerbork Synthesis Radio Telescope) has a resolution of 12 VLA (Very Large Array) 
of 5" — 15", ATCA (Australia Telescope Compact Array) of 15". The IRAM interferometer at Plateau 
de Bure measures CO RCs by means of the CO(l-O) transition at 2.6 mm, at s spatial resolution of 
3". 

Slide 14 

Gran Telescopio Canarias is a 10.4 m diameter reflecting telescope at the Roque de los Muchachos 
Observatory on the island of La Palma. 

In the summit of Hawaii dormant Mauna Kea volcano there are the twin Keck telescopes, the 
optical and the infrared. The telescopes primary mirrors are 10 meters in diameter and are each 
composed of 36 hexagonal segments that work in concert as a single piece of reflective glass. 

The Large Binocular Telescope is located on 3,300 m Mount Graham in Arizona. The LBT is one 
of the world's highest resolution and most technologically advanced optical telescopes; the combined 
aperture of its two mirrors makes it the largest optical telescope in the world. 

On Paranal ESO operates the ESO Very Large Telescope (VLT) with four 8.2-m telescopes. Each 
of them provides one Cassegrain and two Nasmyth focus stations for facility instruments. 

In measuring optical RCs, the class of the telescope is not always decisive, but of course important. 
Telescopes in the 3-4 meter class are sufficient to provide high-quality optical kinematics for local 
objects. The biggest telescopes are necessary to measure the RCs at high redshifts. 

Instead, in measuring dispersion velocities, new science usually requires the best telescopes, also 
for local objects. 

Slides 15-16 

We measure recessional velocities by Doppler shift, and from these, the rotation of the disk. In the 
process, we obtain the sky coordinates of the galaxy kinematical center, its systemic velocity, the degree 
of symmetry, the inclination angle and the circular velocity V{R). Optical high quality kinematics 
can obtain circular velocity at a resolution < Rd/10- 

Slides 17 

During the 70's there was running the idea of a discrepancy between kinematics and photometry, 
already hinted in Freeman (1970) remarks about the NGC 300 RC. A topical moment was when 
Rubin published 20 optical RCs extended to Rppt that were still rising or flattish at the last measured 
point (Rubin et al., 1980). As a guide, we show the profiles of the stellar disk contribution if this was 
the sole contribution to the circular velocity. 

Let us notice, however, that the actual contribution due to the stars, in view of the presence of a 
central bulge and, in some case, of a not exponential disk, may be inside Rppt flatter than how it has 
been drawn (but never rising). 
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4.1 Slide 18 



The decisive measures that detected a dark component around spirals without doubt were obtained 
from the 21-cm RCs. Notice that HI Une require radio telescopes with adequate sensitivity to detect 
the weak signals, and sufficient angular resolution to resolve the details. Only in the course of several 
decades the equipment improved, to arrive in the late 70' at conclusive results. A further improvement 
occurred later by using interferometers, based on combining the signals of two or more telescopes, and, 
by making use of the principle of earth rotation synthesis. Today maps of the integrated HI distribution 
and the radial velocity field are at a resolution of < 10" . 

Bosma (1978, 1981a, b) collected data and compiled HI RCs for 25 galaxies well extended beyond 
the optical radius Ropt, indicated by a red bar on each curve. One obvious case is NGC 5055 shown in 
the slide. It is immediate that the RCs do not decline in a Keplerian way oc R~^/'^ beyond this radius. 

More quantitatively Bosma & van der Kruit (1979) obtained the surface photometry for these 
objects and found, in their outer parts, a clear increase with radius of the local mass-to- light ratios 
'di./dR • light traces the gravitating mass, this quantity is expected to be constant with radius, equal 
to the average mass -to- light ratio of the stellar disk. 



Slide 19 

The radial Tully-Fisher relationship gives a most straightforward evidence for dark matter in spirals. 
The original paper of the standard TF relation is: Tully Sz Fisher (1977). 

At different galactocentric distances, measured in units of the optical size, i?j = i Ropt (i = 
0.2, 1), it has been observed the existence of a family of independent Tully-Fishcr-like relationships, 
Mband = bi + O'i log V{Ri) with Mhand, the magnitude in a specific band, in detail the (R, I) bands. 

This ensemble of relations, we call the Radial Tully-Fisher (Yegorova & Salucci, 2007), constraints 
the mass distribution in these objects. It shows a large systematic variation of the slopes Ui with 
Ri'. they range across the disk between —4 and —8. Moreover the RTF has a very small r.m.s. 
scatter at any radius. This variation, coupled with the smallness of the scatter, rules out the case 
in which the light follows the gravitating mass, shown in the figure as the black line predicting: 
M/ = consti + 10 log V{Ri) for the case of a constant stellar mass-to- light ratio. However, also in 
the case of a magnitude dependence of the latter quantity, the slopes of the RTF could be different 
from the value of 10, but they should be all equal, at every Ri. 

The RTF rules out mass models with no DM or with the same dark-to-luminous-mass fraction 
within i Rgpt in any galaxy. 



Slides 20-21 

Looking to individual RCs (right) we realize that their profiles are a function of the galaxy magnitude, a 
feature that points towards an universality of the spiral kinematics. Persic. Salucci. & Stel (1996) built 
9 coadded RCs from about 700 individual RCs. The whole /-band luminosity —16.3 < M/ < —23.4 
of spirals was divided in 9 luminosity intervals, each of them with ~ 1500 velocity measurements. For 
each interval the corresponding data were coadded in radial bins of size 0.3 -Rd to build the synthetic 
curves VrotiR/Ropt, Mj) out to ~ 1.2Ropt, Ropt = 3.2 Rd- 

The coadded RCs result free from observational errors and non-axisymmetric disturbances present 
in individual RCs see Fig 1 and Bl of Persic. Salucci. & Stel (1996) and show a very small r.m.s. This 
result was confirmed by a similar analysis performed with 3000 RCs by Catinella, Giovanelli, & Haynes 
(2006) (points in the plot). 



Slides 22 

The URC is built from observations out to 6 Rd- There are not reliable kincmatical tracers outside 
this radius. It is possible, however, to extrapolate it out to virial radius (see Salucci et al. (2007)), the 
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bonafide radius of DM halos, (related to the virial mass by Ryir = 26O(M^jr/(lO^^M0))^/^ kpc) by us- 
ing the galaxy virial mass M^r and the virial velocity = G Myir/Rvir- The former can be obtained 
by weak lensing method (Mandelbaum ct al., 2006, 2009) (left plot in the slide), and by correlating 
the galaxy baryonic mass function dN/dMfjar with the theoretical DM halo mass function clN/dMyi^ 
(right plot). Two different recent works that have applied this method well agree with themselves 
and with weak lensing result. (Shankar et al., 2006; Vale &; Ostriker, 2004; Marinoni k, Hudson, 2002; 
Moster et al., 2010; Conroy et al., 2006). Then, at least in first approximation: 



^10 {M^^r/Msf■^ 

(1 + (M^,,/M,)2-2) 



Md = 2.3 X 10^" , ^9 9^ Mq 



with Ms = 3 X 10 that gives V^r in terms of Md, derived by the (known) inner URC. 



Slides 23 

These results lead to the concept of the Universal Rotation Curve (URC) of spirals. Let us define 
X = R/Rd and L the galaxy luminosity. We find that, for any x and L, the Cosmic Variance of 
V{x, L), i.e. the variance of the circular velocity at a same radius x in galaxies of same luminosity L, 
is negligible compared to the large differences of V{x, L) i) in a galaxy of luminosity L when x varies 
and ii) at the same radius x when L varies. 

Therefore, the circular velocity at a given radius, this radius and the galaxy luminosity, lie on a 
smooth surface. 

The URC is shown in the provided movie2 that should be linked to this slide. 



Slides 24 

The individual circular velocities V{r), so as the coaddcd curves are the equilibrium circular velocities 
implied by the galaxy mass distributions. This allow us to decompose them in different mass compo- 
nents contributions. In detail, the gravitational potentials of a spherical stellar bulge, a DM halo, a 
stellar disk, a gaseous disk (^^ot = (t>h + (t>H + <i>disk + (t>Hi lead to: 



with the Poisson equation relating the surface (spatial) densities to the corresponding gravitational 
potentials. In general, we can neglect the contribution of the bulge, but this can play a crucial role 
in certain cases. HI surface density yields V^ag, usually computed numerically. The contribution from 
the stellar exponential disk is 

2 GMd 2T^fX 



where x = t/Rd-, G and B = IqKq — IiKi, a combination of Bessel functions. 

For the DM component we follow the approach of (a quite general) Burkert empirical profile 
(Burkert, 1995) 

p{r) = 



r + ro) (r2 + r^) 



ro is the core radius and po the central density. Vj^^^^ = GMn/r, Mh = Airr'^ p{r)dr . 



„2 



It is sometime used the pseudo- isothermal density p{r) = po -^^xf^. We warn about this often used 

cored distribution (outside the baryonic regions of galaxies where any cored distribution is equivalent) . 
In fact, it implies a divergent halo mass and a outer density profile inconsistent with that of NFW, in 
disagreement with observations. 

Finally, we can adopt the NFW profile emerging out of N-body simulations. 
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The differences among the different halo profiles are shown in the plot. 

The mass model has three free parameters: the disk mass, and two quantities related to the DM 
(the halo central density and core radius for Burkert halos). These are obtained by best fitting the 
data. 

Notice that also for different Hubble type ellipticals, dSphs and for different techniques (weak 
lensing, dynamical modelling of dispersion velocity) the mass modelling always reduces in fitting the 
total gravitating mass profile M(r) by means of a mass model that includes two components: M/,(Pi, r) 
and Mh{Pi, P2,r) and it has, usually, three free parameters Pi (the mass in stars, the halo mass, the 
halo length scale). 

Slides 25 

We start by mass modelling the URC. 

A Burkert halo + Freeman disk perfectly reproduce the coadded RCs at any luminosity. The 
luminous regions R < 2i?£) of the lowest luminosity objects are dark matter dominated, while in those 
of high luminosity objects, the stellar disk is the main mass component. 

In any object, the dark component increases monotonically its importance on the circular velocity 
as R increases. 

Importantly, the structural DM and LM parameters are related among themselves and with lumi- 
nosity. In the 3D figure we sec that spiral galaxies in the 4D space defined by central DM density, core 
radius, luminosity, fraction of DM at Ropt lie on a curve. On a physical side, we realize that smaller 
galaxies are denser and have a higher proportion of dark matter. 

Slides 26-27 

The relationships between halo structural quantities and luminosity must be investigated also via mass 
modelling of individual galaxies (Kormendy Sz Freeman, 2004) that sometime is helped by assuming 
that the stars dominate the very innermost regions. The results confirm the picture derived from the 
URC and point to the quantity poVq which turns out to be constant in galaxies. Cored halos seem 
indispensable and smaller objects have always more proportion of dark matter. 

Slides 28-29 

The investigation of the density distribution of DM around spirals to check whether it complies with 
the 'raw' ACDM predictions needs careful model independent analysis of individual RCs. 

This has been done by using suitable objects (Gentile et al., 2005), and also large galaxy samples, 
including the recent The HI Nearby Galaxy Survey of uniform and high quality data, with a significant 
number of rotation curves in which non-circular motions are small. 

The available rotation curves axe successfully fitted by a cored-halo -|- the stellar/HI disks. For 
DDO 47 the separate contributions are indicated. Instead, they are badly fitted by the NFW (or 
similar) halos -|- stellar/HI disks model, of which for IC 2575 we plot the separate components. The 
blue region indicates the discrepancy between NFW -|- disks models and actual data. 

The results in these galaxies are quite typical. In addition to poor RC fits the latter mass model 
often leads to very large virial halo masses, and stellar mass-to-light ratio that result too low for 
the observed luminosities and colors of the galaxies. It is also proved that tri-axiality and non- 
circular motions cannot explain such CDM/NFW cusp discrepancy with the observed RC data of 
spirals and LSBs (Gentile et al., 2004; de Blok et al., 2008; Kuzio de Naray et al., 2008; Oh et al., 
2008; Spano et al., 2008; Trachternach et al., 2008; Donato ct al., 2009). 

It is however worth recalling that several physical processes that modify the original cosmological 
NFW halo profiles to meet with the observed ones have been proposed. This discussion is beyond 
the aims of this talk, especially because consensus is far from being reached about the net effect of 
(baryonic) physical processes on the dark matter distribution: some physical processes make it close 
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to the observed (cored) (Governato et al., 2010), while others make the NFW cusp even steeper, e.g. 
(Gnedin et al., 2004). 

5 Dark matter in ellipticals 
Slide 32 

In Ellipticals the surface brightness of stellar spheroid, the main baryonic component, follows a Sersic 
(de Vaucouleurs) law: 

I{R) = h d0^[-hn{R/Ref''' - 1] 

where n is the index of Scrsic defining the degree of concentration of the light and is the half 
light radius. For n = 4 we obtain the well-known de Vaucoulers profile By deprojecting the surface 
density I{R) we obtain the luminosity density j(r) and by assuming a radially constant stellar mass 
to light ratio (M/L)^, the spheroid density Psphir)- The central surface brightness /xq is given by: 
jjLQ = 2.5 log /q + const. 

Slide 33 

To derive the gravitational potential in Ellipticals and from this to infer the dark matter distribution is 
significantly more complicated than in spirals. The former can be obtained, always caveat a number of 
crucial assumptions, from dispersion velocities of stars or planetary Nebulae, from the X-ray properties 
of the emitting hot gas halos that lies around them, or from a combination of weak and strong lensing 
data. 

Jeans equation is a crucial link between observations and data. We infer from the observed motions 
the underlying gravitational potential. The knowledge of the distribution of dark and luminous matter 
in these systems of different morphology, age and formation process, is however indispensable. 

Slides 34 

In ellipticals the kinematics is complex, the stars are in gravitational equilibrium by balancing the 
gravitational potential they arc subject to with 3D motions, whose r.m.s exerts a pressure. Obviously, 
from the motions of stars in a galaxy we cannot measure the radial/tangential velocity dispersions, 
directly linked to the mass profile, but recessional velocities or projected and aperture velocity disper- 
sions, quantities that are a much more indirect tracer of the gravitational potential, as shown by the 
formula in the slide, valid for case of isotropic orbits. The kinematics and its analysis becomes more 
complicated when this assumption is relaxed and/or some rotation is present. 

In this case the full 2D kinematics is absolutely required and we must solve (with some assumption) 
the anisotropic Jeans equation, i.e. that involves also the higher moments of the LVDS at any position 
in the galaxy (Cappellari, 2008; Krainovic et al., 2005). 

Slide 35 

Ellipticals are quite compact objects, with respect to spirals. Stellar kinematics of the spheroid traces 
a very inner region of the DM halo, where likely it has small dynamical effect. Moreover, the analysis is 
far from straightforward. This can be seen when we try to reproduce the (observed) velocity dispersion 
profile Uap{R) with a reasonable model made of a dark halo and a luminous mass component. We 
take a galaxy with a spheroid of lO^^M© with the standard Sersic parameter n = 3, a NFW halo with 
c = 7 and a mass 20 times bigger that of the spheroid. In this mock object the DM has a an dynamical 
important role, and it would emerge if its gravitational potential was probed by the rotation curve. 
However, when probed by Gap-, the nature itself of this kinematical hides the dynamical effect of the 
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dark component (Mamon & Lokas, 2005). For a cored halo or with the presence of anisotropy the 
situation would be worse. 

A similar result comes from the work of Tirct ct al. (2010) in which it is shown that an assumed 
typical NFW DM halo around an elliptical, i.e. with M^.^ < 10^"^ Mq, contributes to the observed 
central aperture velocity dispersion (e.g. to aap{l/2 Re)) in a negligible way: CFap^ni^/'^ Ref < 
(100 km/s)^ << (o"ap)^. For cored halos the contribution is even smaller. Only the smallest ellipticals 
and dwarfs spheroidals in which (70^(1/2 Re) < 100 km/s , the stellar velocity dispersion may probe 
the dark component. 



Slide 36 

Important information on the mass distribution can be obtained from the Fundamental Plane. For 
virialized stable objects one expects the balance between the potential and kinetic energies. 

2 GMdyn Mdyn L Mdyn J. jj 

where a is a velocity dispersion, I is a surface luminosity, R is a, scale, and Mdyn/L is the mass-to- 
light ratio. This implies a relationship between the observed velocity dispersion cjg. Re and surface 
luminosity Jg = Le/{2TrRe). i.e. a (Fundamental ) Plane in the space of these 3 quantities. From the 
virial equilibrium, we expect: R^ = CqZ-^q- 

This relationship is found in ellipticals and SO: in the 3D space of {Re,cro,< le >), where the 
latter is the surface brightness inside Re (Jorgensen et al., 1996) we have: 

log Re = 1.24 log ao - 0.82 Log < I >e 

the scatter is very small, 0.07 dex in log Re 

A similar results is found in the large sample of ellipticals that include 6000 objects of the digitalized 
survey SDSS (Hyde & Bernardi, 2009). 

The very small scatter of objects around this plane suggests some degree of universality in their 
inner mass profile. 

The value of the coefficient a is different from that expected by the virial theorem, a = 2 in the 
case in which ellipticals have the same dynamical mass-to-light ratio. This tilt of the FP can explained 
in a dependence of the stellar mass to light ratio (Mgph/L) on luminosity (spheroid mass). In fact, 
FP residuals are found correlated with stellar population characteristics as line- strength indices and 
are anticorrelated with the abundance ratio a/Fe, (Gargiulo et al., 2009). 



Slide 37-39 

Planetary Nebulae (PN) provide us with the elliptical galaxy kinematics out to many effective radii, 
(Mendez et al., 2009). The use of dedicated instruments, like the Planetary Nebula Spectrograph 
(Douglas et al., 2002) has strongly powered the technique, so that distances comparably to those 
reached by extended HI disks in spirals can be probed (Coccato et al., 2009). 

In each galaxy 200 individual PN radial velocity measurements are usually obtained. By binning 
PN radial velocities in elliptical bins, or in strips over the major and minor axes (see for details 
(Napolitano et al., 2001)), reliable 2D kinematical maps have been derived (Coccato et al., 2008). 
These maps led to independent and high quality projected velocity dispersion measurements out to 
(5 — 8)Re, with a typical spatial resolution of 1/3 Re- 

A systematic study of the velocity dispersions in a dozen of objects has revealed that in the 
most luminous objects the velocity dispersion tends to show a regular, monotonic fiat profile, that 
poses no problem in the subsequent mass modelling, other than its uniqueness (Coccato et al., 2009; 
Napolitano et al., 2010). 

In low luminosity objects, instead, the PNs velocity dispersion profile show a large pseudo-and 
ultra-keplerian decline: crp oc i?" with a < —0.5. To analyse these data obviously requires a full 
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modelling of a very likely complex dynamics, that includes also an appropriate treatment of the 
orbital anisotropies of the tracers. 

For objects with "regular" kinematics the best case analysed so far is NGC 4374 (Napolitano et al., 
2010) where a flat dispersion profile has been derived from about 450 PNes and then modeled by means 
of Jeans analysis. 

As result, a variety of halo model profiles including NFW, some of its variants and the cored URC 
halo, in cooperation with a spheroid with typical mass-to-light ratios, well fit the data, confirming the 
presence a very massive dark halo. Notice that in this elliptical the kinetical energy oc V^^g has both 
rotational and random motions component: Vrms = O^rot + '^p)^^'^ ■ 

Slides 40-42 

Gravity bends light and this can be used to weigh the galaxy mass, e.g. the pioneering work by 
(Schneider, 1996). We know that a background round galaxy (of semi axis a/b = 1) whose line of 
the sight passes at distance Rp from a the center of a galaxy put at a convenient distance (half way 
between us the background) will be seen with a (very small) ellipticity a/b—1 ~ 10^^. The ellipticity 
of a galaxy image is an unbiased estimate of the local shear. When averaged on a very large number of 
source galaxies whose lines of sight all pass around a lens, we can meaningfully measure the tangential 
shear jt of the lens. The lens equation relates jt with the distribution of matter in the lensing galaxy: 

7i = (E - S(i2))/S„ 

where S(i?) = 2 piR, z)dz is the projected mass density of the object distorting the galaxy image, 

at projected radius R and S(i?) = xT,{x)dx is the mean projected mass density interior to the 

radius R. = ^^ ^^^ , where Dg and Di are the distances from the observer to the source and lens, 
respectively, and Dis is the source-lens distance. The above relations directly relate observed signals 
with the underlying DM halo density. A very detailed explanation of the weak lensing phenomenon 
can be found in the second chapter of Schneider, P., Kochanek, C. S., &: Wambsganss, J. (2006). 

Let us stress that we discuss this effect in the section of elliptical in that it has been applied mostly 
for these objects, however there are applications (even in this Review) for Spirals. 

The DM distribution is obtained by fitting the observed shear with a chosen halo density profile 
with 2 free parameters. The method is clean because it is free from the baryonic contribution, but 
the signal is intrinsically weak. Mandelbaum et al. (2009, 2006) measured the shear around 1.7 x 10^ 
isolated galaxies of different luminosities and Hubble Type, by using 3 x 10^ SDSS galaxies as sources. 
The measurements extended out to 500 — 1000 kpc reaching out the virial radius of the lens galaxies. 
Although the structural DM parameters were determined with quite large uncertainties, the NFW 
and Burkert halo profiles, not much different at outer radii, both agree with the data. No relevant 
difference in the DM distribution around ellipticals or spirals emerged. 

Very importantly, the virial halo masses so determined are found to correlate with the galaxy 
luminosities. When this relation is coupled with the well known mass vs. light relationship of the 
stellar component, it leads to the fact that the amount of baryons M^aj. today in a galaxy is a function 
of M^ir , the mass of the halo surrounding the galaxy. We consider this as the main relationship of the 
galaxy formation process. It is interesting to note that M^ir/M{^ar >> 7, the cosmological value. 

Noticeably, when we fit the tangential shear with a Burkert profile, we obtain the same val- 
ues of halo structural parameters po , ro derived by mass modelling the URC or individual galaxies 
(Donato et al., 2009). 

Slide 43-44 

Koopmans et al. (2006) performed a joint strong lensing -|- stellar-dynamical analysis of a sample of 
15 massive early-type galaxies selected from the SLAGS Survey. Given the smallness of the radial 
range investigated, the total mass density was assumed a power law without loss of generality. The 



13 



spheroid component, was represented for simplicity, as an Hernquist profile, very similar inside Rg to 
a Sersic one. Velocity dispersions at Re provided the additional constraint to the mass distribution, 
that was investigated through Jeans Equation. 

They found, independently of the DM density profile in the region under investigation, that these 
massive early-type galaxies have remarkably the same inner total density profiles {ptot = Ph + Psph «^ 
r~^). The figure shows the logarithmic density slope of SLACS lens galaxies as a function of (normal- 
ized) Einstein radius. Moreover, they found that the stellar spheroid accounts for most of the total 
mass inside Rf,. 

In both cases, given the smallness of the dark component inside R^, the small radial range, the 
not negligible observational error and the relevance of the assumptions taken, it is not possible to 
constrain the actual profile of the dark halo. 



Isolated Ellipticals have a X-ray emitting halo of regular morphology, that, extends out to very large 
radii, providing us with a reliable mass distribution. The gravitating mass inside a radius r, M(r) 
can be derived from their X-ray flux by assTiming that the emitting gas is in hydrostatic equilibrium. 
The observed X-ray surface brightness, gives, by means of the well known j5-model, the gas density 
distribution. From the hot gas density and temperature profiles (Fabricant, Rybicki, & Gorenstein, 
1984), we have: 



where Tg is the gas temperature at radius r, pg is the gas density, k is the Boltzmann's constant, G is 
the gravitational constant, ji is the mean molecular weight ji = 0.62 (Ettori &; Fabian, 2006), and rup 
is the mass of the proton. 

Nagino & Matsushita (2009) analyzed a sample of 22 objects with XMM-Newton and Chandra 
data out to 10 Re- They confirmed the presence of DM in any object, but this component was found 
dynamically important only outside Re, where the mass-to- light ratio is found to start to increase. 
Furthermore, there is some hint that, outside this radius, the halo mass increases with radii quite 
steeply suggesting a cored DM halo density, for which the dark mass increase rate between and r^. 



The gravitating matter in form of stars is a crucial quantity. It defines the efficiency of the process 
of forming this (luminous) component from the cosmological hydrogen. Usually, it is obtained as result 
of the kinematical mass modelling, but in some case its resulting value has a very large uncertainty. 
This, in turn complicate the full investigation on the Dark Matter distribution n. Moreover, in low 
luminosity and DM dominated objects, the kinematics cannot provide a reliable estimate of the stellar 
disk mass. 

It is well known that the galaxy luminosity is related to its stellar content. Hence the straight- 
forward approach to obtain the stellar mass of a galaxy is to model its luminosity in terms of age, 
metallicity, initial mass function. This modelling is provided by Stellar Population Synthesis models. 
This method is presented in this section since it has important applications in Ellipticals, but it is 
applied also to Spirals and dSphs. 

The observed Spectral Energy Distribution (SED) of a galaxy, or selected colour indices or ab- 
sorption lines, is fitted with models calculated under different assumptions regarding the IMF, the 
star formation history, the metallicity, etc. The same modelling naturally provides the amount of 
stellar mass that is implied by the different fits (for a classical paper see Tinsley (1975). In practice, 
degeneracies between age, mctallicitics and dust content make it sometimes hard to select the actual 
physical model for the galaxy light just in terms of fit goodness. Observations are decisive to solve 
these intrinsic uncertainties. 



Slide 45 




Slide 46 
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While the best approach is obviously to fit the whole spectral energy distribution, as this allows 
a better handle of elusive degeneracies, often not all observational bands are available for galaxies. A 
simplified approach was suggested by Bell &: dc Jong (2001) who noted that rather simple relationships 
exist between mass-to-light ratios in certain bands and some colour indices of population synthesis 
models (from various codes). These relations have been calibrated using spiral galaxies. In the plot 
there are two different M/L's as a function of (B-R), for a sequence of exponentially declining star 
formation rate models of age 12 Gyr obtained by using a variety of Stellar Population Synthesis models. 
The red end of the lines represent a short burst of star formation, the blue end represents a constant 
star formation rate model. The different models used are: Bruzual &: Chariot, Kodama &; Arimoto, 
Schulz et al. and PEGASE by Fioc Sz Rocca-Volmerang all with a Salpeter Initial Mass Function. 

Recent work considers in great depth the TP-AGB phase of stellar evolution while building the 
galaxy SED (Maraston, 1998, 2005) which were shown to provide a better fit to the spectral energy 
distribution of high-redshift galaxies (Maraston et al., 2006), allowing to better determine their stellar 
masses. The resulting mass-to-light ratios, as a function of galaxy color are shown in the slide. 

Slide 47 

The photometric estimates of stellar disk/spheroid masses are tested with respect to the kinematical 
estimate of the mass in stars. 

In Salucci et al. (2007) the disk masses of 18 spiral galaxies of different luminosity and Hubble 
Type have been derived both by mass modelling their rotation curves and by fitting their SED with 
spectro-photometric models. The two estimates result in perfect proportionality and the agreement 
between the two different estimates is less than 0.15 dex. 

In Elliptical galaxies the determination of dynamical mass inside Rg can be performed by as- 
suming that a De Vaucoulers spheroid is the main component then: M^y^ = 1.5 x 10^'^((Te/200 
km /s)'^ {Re/ kpc)M_odot This dynamical estimate of the "stellar" mass is found to well correlate with 
the red luminosity, i.e. with a measure of the photometric mass estimate (Shankar & Bernardi, 2009). 

Grillo ct al. (2009) studied a sample of Ellipticals with Einstein rings and derived the total pro- 
jected mass enclosed within them. Then, by using the SDSS multicolor photometry fitted the spec- 
tral energy distributions (SEDs) (ugriz magnitudes) with a three-parameter grid (age, star-formation 
timescale, and photometric mass) of Bruzual & Chariot's and Maraston's composite stellar-population 
models, and various initial mass functions (IMFs). The two different mass estimates agreed with 0.2 
dex. 

A similar result has been obtained by studying a sample of ellipticals with detailed axisymmetric 
dynamical modelling of 2D kinematics (Cappellari et al., 2006). 

Caution: the disagreement among dynamical and spectro-photometric estimates is about 0.15 dex, 
this value is small for cosniological implications, i.e. to build a stellar mass function from a luminosity 
function, but it is extremely large for mass modelling. The induced error, if we assume {M^/LY^^ in 
estimating the DM profile 

PH{r) = ^d{Mayn{r) - {M,/ Lf^ L{r))/dr 
can totally err the derived density profile. 

6 Dark matter in dwarf spheroidals 
Slide 50 

Dwarf spheroidal (dSph) galaxies are the smallest and least luminous galaxies in the Universe and so 
provide unique hints on the nature of DM and of the galaxy formation process. Early observations 
(e.g. Aaronson (1983)) of stellar kinematics in the brightest [Ly ~ lO^^^Ly^©) of the Milky Way's 
dSph satellites found that they have central velocity dispersions of ~ 10 km s~^, larger than expected 
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for self-gravitating, equilibrium stellar populations with scale radii of ~ 100 pc. The recent discoveries 
of ultra-faint Milky Way satellites (e.g. Willman et al. (2005); Belokurov et al. (2007)) extend the 
range of dSph structural parameters by an order of magnitude in (luminous) scale radius and by three 
orders of magnitude in luminosity. 

The distribution of total luminosity and of the central stellar velocity dispersion in star clusters 
and in dSphs overlap. The faintest galaxies have approximately the same values of these physical 
parameters of star clusters, with galaxy luminosity extending down to IO^Lq. In dSphs the half- 
light radii, however, are significantly larger (hundreds of parsccs) than those of star clusters (tens of 
parsecs). This leads, through the virial theorem, to significantly larger masses in dSph. 

The values of central and global mass-to-light ratios for the gas-poor, low-luminosity, low surface 
brightness satellites of the Milky Way are high, up to several hundred in solar units, making these 
systems the most dark matter dominated galaxies in the local Universe (Mateo, 1998). 



DSphs arc old and apparently in dynamical equilibrium, with no dynamically-significant gas. They 
contain a (small) number of stars, which provide us with collisionless tracer particles. Aaronson (1983) 
measured the velocity dispersion profile of Draco based on observations of three carbon stars, finding a 
mass-to-liglit ratio an order of magnitude greater than that found for galactic globulars. Mateo (1997) 
provided the first dispersion velocity profile of Fornax. With the new generation of wide-field multi- 
object spectrographs on 4-8m telescopes, it is now viable to obtain sufficient high-quality kinematic 
data to determine the gravitational potentials in all the Galactic satellites. 



Walker et al. (2007) published high quality velocity dispersion profiles for seven Milky Way satellites. 
They are generally fiat out to large radii, in terms of the half light radii that determine the sizes of 
the baryonic matter in these objects. These data indicate that almost all of their gravitating mass is 
unseen, i.e. in a DM halo. 



The mass modelling of these object is not immediate, in that the distribution of the stars, that are the 
tracers of the gravitational fields, must be known with precision. Moreover, the relation between mass 
and dispersion velocity is not a straightforward one: in a collisionless equilibrium system that of the 
stellar spheroids in dSph, the Jeans equation relates the kinematics of the tracer stellar population and 
the underlying (stellar plus dark) mass distribution. Assuming spherical symmetry, the mass profile 
can be derived as 



where ar{r), (3{r) and v{r) are, respectively, the stellar velocity dispersion component radially toward 
the center of the mass distribution, the velocity anisotropy and the stellar density distribution. 

Disadvantages in modeling these systems are: i) it is required to take some assumption on the 
velocity anisotropy and ii) derive the stellar density profile with precision, in order to fit the observed 
velocity dispersion profile and obtain a reliable DM distribution. 

An advantage is that their luminosities are so small that we can always exclude that baryons play 
a (complicated) role in the mass modelling. 

The surface brightness profiles are typically fit by a Plummer distribution (Plummer, 1915) 



where R is the projected radius, vq is the central surface density and is stellar length-scale Notice 
that, if we want to derive the DM density profile (not its amount) we must know with good accuracy 
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Slide 53-54 
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the distribution of the stellar component, since this quantity enters in the determination of the mass 
profile in a non- linear way. In Elliptical and Spirals, it is insread extremely rare that, given a set of 
kinematical data, two different stellar distributions but both compatible with the measured surface 
photometry, yield, in the mass modelling, different DM profiles. 

The analysis of the mass distributions obtained from isotropic Jeans equation analysis of six dSphs 
by Gilmore et al. (2007b) favor a cored dark matter density distribution. 

However, ? applying the Jeans equation to eight brightest dSphs and relaxing the isotropy as- 
sumption found that the dispersion velocities of these objects can be fit by a wide range of different 
halo models and velocity anisotropics. In particular cusped and cored dark matter mass models can 
fit the dispersion profiles equally well. 

Slide 55 

A resolutive fact could be that dSphs cored distribution structural parameters surprisingly agree with 
those of Spirals and Ellipticals: Donato et al. (2009) have shown that a Burkert halo density profile 
reproduces the dispersion velocity of six Milky Way dSphs and that the derived halo central densities 
and core radii are consistent with the extrapolation of the relationship between these quantities found 
in Spirals. This intriguing result, beside supporting cored distribution in dSph, could point to a 
common physical process responsible for the formation of cores in galactic halos of all sizes or to a 
strong coupling between the DM and luminous matter. 

Slide 56 

The dark+ stellar mass within the effective radii of dSph galaxies have been suspected for some years 
of showing a remarkably small range (Mateo, 1998). The available data seem to be consistent with an 
(apparent) minimum dark halo mass, supported also by the presence of a strong relationship between 
M^/Lv and L (Wilkinson et al., 2006; Gilmore et al., 2007a). 

More explicitly, a relationship is found for which, within 300 pc, all object seem to possess the 
same mass, of the order of WMq, (Cote et al., 1999; Strigari et al., 2008; Gilmore et al., 2007b). 
This relationship maybe be not intrinsic, but due to the universality of the profile of the density of 
the halos around galaxies. 

On the other hand, if all dSphs are embedded within an "unique" dark matter halo, this could 
explain the finding that their masses increase as power-law in all objects, M(r) oc r^'^ (Walker et al., 
2009). 

Slide 58 

DM halos around galaxies form a family of dark spheres with their virial mass being the order pa- 
rameter. It is possible to work out an unified scheme for their structural characteristics, that, let us 

remind, are derived by means of several different methods in galaxies of different luminosities and 
morphologies. From all available data a picture emerges: DM halos around spirals, LSBs and dSphs 
and Ellipticals can be well represented by a Burkert profile with a core radius ro and a central density 
Po- 

Very intriguingly: the halo central surface density oc po ''o is found nearly constant and inde- 
pendent of galaxy luminosity. This result come from the mass models of 1) about 1000 coadded 
rotation curve of Spirals, 2) individual dwarf irregulars and spiral galaxies of late and early types 3) 
galaxy-galaxy weak lensing signals 4) kinematics of Local Group dwarf spheroidals. 

It is found: 

log{ro Po) = 2.15 ± 0.2 

in units of log{MQ/pcP). 

This results are obtained for galactic systems spanning over 14 magnitudes, belonging to different 
Hubble Types, and whose mass profiles have been determined by several independent methods. In the 
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same objects, the constancy of po vq is in sharp contrast to the systematical variations, by 3-5 orders 
of magnitude, of galaxy properties, including po and central stellar surface density. 

Slide 59 

A connection between the a-priori independent distributions of dark and luminous matter in galaxies 
emerges strongly. The velocity (mass) profile, despite being the sum in quadrature of both components, 
can surprisingly be fully represented by a function that includes only one of them. In fact, for Spirals, 
(but likely in galaxies of other Hubble Types) we have : V{r) = Fi,ar{r / Rd, M^) (M/ is the I- band 
magnitude, left- side plot), but, at the same time, we also have: V{r) = Fmr{r/Rvir, Myir) (right-side 
plot) (Salucci et al., 2007). 

In this last case the baryonic component breaks the mass invariancc of the RC profile due to the 
DM component, but in a mass dependent way. The differences of the inner RCs profile of high and 
low mass objects are due to differences in the fractional content of baryonic matter, but mysteriously, 
they, in turn, depend on the halo mass. 

It may be worth recalling that similar result (in qualitatively way) are present in the predictions 
of ACDM, where NFW halos coupled with a stellar disks formed according to the Mo et al. (1998) 
theory leads to RCs Universal profile, function of the halo virial mass. 

Slide 60 

Moster et al. (2010) used a statistical approach to determine the relationship between the stellar 
masses of galaxies and the masses of the dark matter halos in which they reside. They obtained 
(following Shankar et al. (2006)) a stellar-to-halo mass (SHM) relation by populating halos and sub- 
halos in an N-body simulation with galaxies and requiring that the observed stellar mass function be 
reproduced. The existence of this non linear relation is not trivial: M^, the mass in stars, does not 
scale proportionally to Myir and it is much smaller than that expected by the Cosmological value (blu 

line) i^baryons/^DM- 

Walker et al. (2010b), subtracted in spirals and in dSph the baryonic contributions to the overall 
gravitational potential to derive the contribution to rotation curves due to dark matter halos at 
intermediate radii. They found that, in physical units, these RCs are remarkably similar, with a mean 
curve given by 

(fc^ = ^-^^-0:^9 + 0.2 log (r/kpc) 

Evaluation at specific radii of this Universal Profile immediately generates two results discussed 
above: a common mass for MW dSphs at fixed radius and a constant DM central surface density for 
galaxies ranging from MW dSphs to spirals. This result is very remarkable: in the Universal Functions 
discussed above as the URC and the FP, the radial unit is always measured in terms of the stellar 
galaxy size and not, as here, in physical units. 

However, in dSphs it is possible to estimate "the circular velocity" only at one radius in each 
galaxy (Walker et al., 2009), namely the half light radius, whose definition itself,however, may be not 
trivial (Walker et al., 2010a). 

This very recent result needs more investigation, though it seems coherent with the phenomenology 
of the DM distribution around galaxies. 

Slide 61 

Let us define Re as the effective radius encompassing halfoi the mass of stellar component and provides 
so a the measure of the size of the latter. At Re in dispersion-supported galaxies, within a large range 
of luminosities, the dynamical mass-to-light (M/L) ratios can be obtained (Wolf et al., 2010). At the 
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high mass end, these estimates come from the FP, while at the low mass end, they come from Jeans 
modelling of dwarf spheroidal galaxies. 

Confronting the emerging mass-to-light ratios, galaxies follow a characteristic U-shaped curve 
(Wolf et al., 2010), where M/L is minimized for ~ L* galaxies and rises at both low and high masses. 
This can be seen in the plot, where the /-band M/L ratio within the half-light radius is plotted against 
the mass contained within the half-light radius. 

One can interpret this as the inverse efficiency with which galaxies retain and turn their initial 
quotient of baryons into stars : 1/7 M^r- There must be two mechanisms that prevent galaxies from 
turning their initial baryons into stars, one of which operates most efficiently at low masses, switching 
on strongly below M(i?e) < 10^ Mq, and one of which operates most efficiently at high mass, gradually 
becoming more effective at masses above M(i?e) > 10^^ Mq. Some suggested mechanisms for the 
former include heating from re-ionization, supernova feedback, and tidal and ram pressure stripping 
for those galaxies that are satellites, while proposals for the latter mechanism are feedback from active 
galactic nuclei, virial heating of gas in the deep potential well, and supernova feedback. 

It is interesting that, globular clusters, whose masses can be calculated using the same Jeans 
modelling methods as those of dwarf spheroidal galaxies, are unambiguously not part of this relation. 
They have roughly constant M/L ratios, consistent with that expected from their stellar populations. 
Systems containing dark matter and those without dark matter are therefore easily distinguished. 

It is worth recalling that the same U shape is seen in disk systems at Ryir (Shankar et al., 2006) 
and it is shown as yellow region in the plot. 

7 Dark matter searches 
Slide 62 

To introduce the various candidates proposed as the Dark Particle and to discuss the present status of 
their possible detection would need a complete (long) Presentation Review. As matter of fact we deal 
these subjects in a brief section of our Review. This is certainly insufficient however this Presentation 
is likely to be for some the first (serious) contact with the DM world, and therefore we must consider 
here such issue. 

Slides 63-65 

DM is actively searched in a direct and indirect way. The distribution of dark matter in our galaxy, 
so as in nearby ones, plays an important role in it. 

Among the many candidates suggested for dark matter but as yet undetected, one among the most 
compelling possibility is that the dark matter is comprised of Weakly Interacting Massive Particles, 
or WIMPs, a general class of particles, thermally produced in the hot early universe, which drop out 
of equilibrium when they are not relativistic. Bertone et al. (2005) discuss the theoretical motivations 
and introduce a wide array of candidates for particle dark matter. 

Its density today is inversely proportional to their annihilation rate, and in order to represent the " 
CMB observed" 25% of the critical density of the Universe, the particle's annihilation cross section 
should be typical of electroweak-scale interactions, hinting at physics beyond the Standard Model. The 
existence of such a particle has plenty of theoretical justification. In fact, several extensions to the 
Standard Model lead to WIMP candidates. One of them is Supersymmetry (SUSY), which extends 
the Standard Model to include a new set of particles and interactions that solves the gauge hierarchy 
problem, leads to a unification of the coupling constants, and it is required by string theory. The 
lightest neutral SUSY particle, the neutralino, is thought to be stable and is a natural dark matter 
candidate. 

Direct detection of WIMPs requires sensing scattering between WIMPs in the halo and atomic 
nuclei in a terrestrial detector. Despite the expected weak-scale coupling strength for such interactions, 
terrestrial experiments are made feasible by the fact that: i) the local DM mass-energy density is 10^ 
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times higher than that of the universal average. In detail, let us stress that ph{Rp)), the DM density 
at the Sun location, has been recently estimated it in an independent way by Salucci et al. (2010): 
Pdm{Rq) = (0.43 lb 0.1) Gev cmT^ a value somewhat larger than that of 0.3 GeV cm~^ usually 
assumed, ii) coherence effects may amplify the interaction rate (e.g. spin-independent nuclear elastic 
scattering is proportional to A'^ the atomic weight of the target). 

Even though direct detection may be possible, it is by no means easy. The energy of the resulting 
nuclear recoils is likely of order 20 keV, where electromagnetic backgrounds dominate by many orders of 
magnitude. Nuclear collisions due to neutrons also pose a serious obstacle, as they may exactly mimic 
those due to WIMPs. Hence, the fundamental challenge of a direct search is to identify extremely 
rare DM-induced nuclear recoil events from a vast array of spurious signals. This can be achieved by 
combining low-radioactivity materials and environments with robust background rejection of electron 
recoil events. These detectors must also be situated in deep underground laboratories to shield from 
cosmic-ray-induced backgrounds. 

The dark matter community has developed a broad range of techniques to address these challenges. 
Many different experiments are currently pursuing one or more of the multiple possible detection 
channels. The most common technique is to look for ionization and/or scintillation of the detector 
medium (e.g. Nal, Csl, Ge, liquid/gaseous Xe or Ar) resulting from a WIMP-nucleon collision. Other 
experiments look for phonons, bubble nucleation, directionality effects. 

About 50% of all direct search experiments use noble liquids (including XENON, LUX, XMASS, 
ZEPLIN, WARP, ArDM, and DEEP/CLEAN). Noble hquids, especially Xe and Ar, are promising as 
a detection medium as they are efficient scintillators, are chemically stable, and have relatively high 
stopping power due to high mass and density. 

At present, the upper limit on the WIMP-nucleon cross section approaches 10~^^cm^, which is 
well into the region of parameter space where SUSY particles could account for the dark matter. The 
next 2-3 orders of magnitude represent a particularly rich region of electroweak-scale physics. The 
next generation of direct search experiments are poised to verify or reject a large portion of the current 
theoretical explanations of WIMP-like DM. 



Slides 66-68 

DM particles makes the galactic halo where they would continuously annihilate into quark or gauge 
boson pairs leading eventually to rare antimatter particles and high-energy photons and neutrinos. The 
indirect serches of DM are based on astrophysical observations of the products of DM self annihilation 
or decay. Given the known long lifetime of the DM, the signal for decay products is suppressed 
for heavy candidates (due to the combination of low number densities and long lifetime) leaving in 
most of the cases only the self-annihilation as most sensitive possible source of a signal. Antimatter 
particles such as antiprotons and positrons are interesting because they are quite rare in Nature. Any 
(additional) production through DM annihilation will induce an excess in the energy distributions of 
the flux observed at the earth. However, let us stress that conventional astrophysical mechanisms in 
the galaxy also produce antiprotons and positrons. 

The indirect searches of DM are based on astrophysical observations of the products of DM self 
annihilation or decay. Given the known long lifetime of the DM, the signal for decay products is 
suppressed for heavy candidates (due to the combination of low number densities and long lifetime) 
leaving in most of the cases only the self-annihilation as most sensitive possible source of a signal. 

In the case of searches via gamma ray observation, the expected flux in a detector on Earth is 
given by: 

where is the photon energy, is the DM particle mass, Aip is the detector opening angle, {p'v)a,nn 
is the mean annihilation cross section times the relative velocity (of order 10~^^ cm^ s~^ for cold WIMP 

relics from abundances constrains), Bf indicates the branching fraction in a given channel /, is 
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the photon spectrum for a given annihilation channel which depends on the DM model and can have 
both continuum and discrete lines contributions, p is the DM density and the integrals are along the 
line of sight and over the detector angle. 

The quadratic dependence on p suggest that the preferred targets for indirect searches are the 
places with higher DM concentrations, like the centre of galaxies, satellites or galaxy clusters. It must 
be noticed however that the galactic centres are very often sources of strong activities due for example 
of the presence of black holes or other compact objects enhancing the overall background. Moreover the 
large uncertainty on the DM density profile directly reflects on the flux predictions making the searches 
extremely difficult (although possible enhancements due to local DM over-densities are possible). 

DM particles can also generate high-energy photons. The Fermi satellite for example maybe able to 
indirectly detect DM particles looking at high energy photons.As example let us consider DM species 
at 40 GeV WIMP from a supersymmetric theory. That particle is supposed to annihilate into bottom- 
antibottom quark pairs. In this scenario, we show in the slide the flux map generated after five years 
of satellite operation. The map features the DM signal to noise ratio. The DM signal is the number 
of gamma rays produced towards each pixel by the DM distribution. The smooth DM distribution 
and the heaviest and closest DM substructures (bright spots on the maps) , as expected in AC DM are 
considered. 

The noise includes the galactic diffuse emission mostly arising from the interactions of cosmic rays 
protons and helium nuclei on the HI and H2 gas components of the Milky Way. The two maps refer to 
different set of simulations: Aquarius (left panel) and Via Lactea II (right panel). It is then obvious 
how the DM density distribution in galaxies pH plays a role, in particular, according whether in the 
innermost Ryir/20 region it is cuspy pH{r) oc or constant. The predicted flux can differ by several 
orders of magnitude. Furthermore, one should also take into account the flux enhancements due to 
local DM over-densities i.e. sub-halo, but this is a very open question. 

Another possibility is to pursue indirect search by looking at charged particles such as positrons 
or antiprotons, in this case however, the galactic magnetic field is such that the direction of arrival of 
the particle does not reflect the production point and the only observable is an excess of antimatter 
with respect to the expected background due to ordinary cosmic rays (which also suffer from big 
uncertainties). 

Among the most important facilities for indirect signals we have: XMM Newton and Chandra, 
Integral, Compton Gamma Ray Observatory, AGILE, Fermi-LAT Observatory, CANGAROO, HESS, 
MAGIC, VERITAS, AMANDA, ICECUBE, ANTARES, PAMELA and AMS. 

Up today several claims of indirect DM detection has been made, sometimes in conflict with each 
other and all of them competing with other astrophysical explanations. The most signiflcant are: the 
positron excess measured by HEAT and PAMELA, the 511 keV line excess measured by INTEGRAL, 
the EGRET Diffuse Galactic Spectrum and the the so called WMAP Haze (an excess of microwave 
emission around the centre of the Milky Way) see Morselli (2010) for the current status on these 
searches. . 

8 What we know? 
Slides 69 

The distribution of DM halos around galaxies shows a striking and complex phenomenology. We 
believe that it is giving invaluable information on the Nature itself of dark matter and on the galaxy 
formation processes. 

Semi-analytical models as Baryon -)-DM simulations should reproduce the empirical scenario: a 
shallow DM inner distribution, strong relationships between the halo mass and 1) central halo density, 
2) baryonic mass, 3) half-mass baryonic radius and 4) baryonic central surface density, a constant 
central halo surface density. In any case, the mass discrepancy in the World-Universe of galaxies is a 
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clear function of radius, total baryonic mass and Hubble Type 

Mb ~ l + 72(M6,T)r2 
that very unlikely may arise as effect of a change of the gravitation law. 
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DARK MATTER IN GALAXIES 

NAME 
INSTITUTE 



Outline of the Review 

Dark Matter is main protagonist in the Universe 




This review focus: Dark Matter In Galaxies 



The concept of Dark Matter in virialized objects 
DarkMatterin Spirals, Ellipticals dSphs 
Dark and Luminous Matter in galaxies. Gliibal properties. 
Phenomenology af the mass distribution in Galaxies. 
Implications for Directand Indirect Searches 




The Realm of Galaxies 



The range of galaxies in magnitudes, t^pes and central surface 
densities : 15 mag. 4 types, 1S mag arsec^^ 



Central surface brightness vs galaxy magnitude 



What is Dark Matter ? 

In a galaxy, the radial profile of the gravitating matter IVI(r) does not 
match that of the luminous component MlO'). 

A MASS VE DARK COMPOrjEr-TT le then inlrodiced to account for the disagreement 
Its prof Je M^(r) must obey: 

dtoeA/(r) ML(r)dhe.MLlr) Majr) dhgMyjr) 
<ik)gr Mir) ilogr M(r) dloer 



M (r), M, (r), dlog M. [r)/dlog r observeiJ 



THEORY AND SIMULATIONS 



The DM phenomenon ci 



iwestiysted cinly rl w 




Spirals -. stellar dlsK +b jIqs ■^HI disk 
Elllptlcals&dwarfs E: stellsr spheroid 



/\CDM Dark Matter Density Profiles from M-bady simulations ' 
The density of virialized DM halos of any mass is empirically described at all 
timas by an Universal profile (Nauflrrn + ^i:, 97, NFW). 



Aquarius N-Body simulaflons, highest mass resolLition to date. 

Density distribution: the EInasto Law indistinguishable by NFW 



lnl@w lhiiil®s if®™ siirodl ®¥@D¥® 



eailler have larger ouerden&itiai 
Today mean halo density inside 
R„,. = lOQ 0. 






SPIRALS 




Circular velocities from spectroscopy 

-Optical emisaion lines (Ho- Na) 

- Neutral hydrogen {HI>carbon monoxide (CO] 






Stellar Disks 






IVI33 disk very smoolfi, 
IruriuatBCl at 4 :jt;a|y-lerig[lw 


□ 




B 


: "'Us- 


NCC3D0 sKponenlial disk 
forat least 10 itaie-ieiigtiii 



Gas surface densities 



i^D ierght scale of 1lie ( 



Flatfish ra 
Deficiency In the centre 
Evtenrled to (8 - 401 Rp 



Follo'A3 the stellar di 



Wonga Blitz (2002) 




ROTATION CURVES 



1 








iZ 




i ™ 


f 



Figure 1: Zoomed out of the first 
shown here to give a ghmpse of it. 



15 shdes of the Presentation Review, 
These shdes can be downloaded at 



http : //www. sissa. it . ap/dmg/dmaw_presentation.html 
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